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I. INTRODUCTION 

The coupling between a quantum system and its environ- 
ment typically causes decoherence, which is detrimental in 
quantum information processing (QIP) as it results in com- 
putational errors [1|. Over the years, various ways of sup- 
pressing quantum decoherence have been explored; see, e.g., 
Ref. Q. The methodology we study here is dynamical de- 
coupling (DD), which utilizes sequences of strong pulses to 
decouple the system from its environment . Recently, an 
optimal DD pulse sequence was discovered for the suppres- 
sion of pure dephasing or longitudinal relaxation of a qubit 
coupled to a bath with a hard high frequency cutoff: Uhrig 
DD (UDD) I8UI31 In UDD, the instants tj ( j € {1,2... N}) 
at which 7Y instantaneous it pulses are applied are given by 
tj = TXj, where T is the total time of the sequence, and 



By optimal it is meant that each additional pulse suppresses 
dephasing or longitudinal relaxation to one additional order in 
an expansion in powers of T, i.e., TV pulses reduce dephas- 
ing or longitudinal relaxation to 0(T N+1 ). Rigorous perfor- 
mance bounds were established in Ref. [14|. In this work we 
derive performance bounds for more general pulse sequences. 

A near-optimal way to suppress general single-qubit deco- 
herence, as opposed to only pure dephasing or pure longitudi- 
nal relaxation, is the quadratic DD (QDD) sequence lfT5l4T9l . 
A QDD sequence is obtained by nesting two UDD sequences 
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of pulses which are orthogonal in spin space. When these two 
UDD sequences comprise the same number of pulses, 7Y, a 
QDD sequence of (N+ 1) 2 pulses will suppress general qubit 
decoherence to 0(T N+1 ), which is known from brute-force 
symbolic algebra solutions for small N to be near-optimal 
JT3J . It has been analytically proven that when the two UDD 
sequences comprise different numbers of pulses, Ni and N2, 
a QDD sequence of (Ni + 1)(A?2 + 1) pulses will suppress 
general qubit decoherence to 0(T min( - Nl - N ^ +1 ), and this is 
universal, i.e., holds for arbitrary baths and system-bath inter- 
actions JT61 [T8l [191 - Moreover, the dependence of the sup- 
pression order of each single qubit error type (dephasing, bit 
flip, or both) on Ni and N2 has also been established in detail 
both in terms of analytical bounds 1 18 1 and numerical simula- 
tions (1711201 . 

When the nesting process of UDD sequences is continued, 
one has the nested UDD (NUDD) sequence, which can pro- 
tect multiple qubits, or general multi-level quantum systems, 
against general decoherence |[T6l [191 . It has been proven that 
NUDD is universal, and will suppress general multi-qubit de- 
coherence to 0(T waini ^ N '^ +1 ), where Ni are the orders of the 
UDD sequences being nested [19|. However, it is known that 
NUDD is sub-optimal ifTrSl . 

Our main results in this paper are analytical upper bounds, 
for QDD and NUDD, on the trace-norm distance between the 
states of the DD-protected qubit or qubits, and the unperturbed 
qubit or qubits, given as a function of the total evolution time 
and the norms of the bath operators. Under the assumption 
that the bath operators have finite norms, the upper bound for 
NUDD shows that the trace-norm distance can be made arbi- 
trarily small as a function of the minimal decoupling order N 
of the UDD sequences comprising the NUDD sequence. In 
the QDD case, a tighter bound is obtained by having different 
decoupling orders for different types of decoherence errors, 
using the results of Ref. (T8l . 
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The structure of this paper is as follows. We develop our re- 
sults for QDD and NUDD in parallel, always starting with the 
simpler case of QDD. We first review the QDD and NUDD se- 
quences in Section [TT] We also derive bounding series for the 
two sequences in this section. In Section III we use the bound- 
ing series in order to find explicit upper bounds on the dif- 
ferent single-axis errors for QDD, and similar explicit upper 
bounds on different error types for NUDD. These results are 
then used in Section [IV] to derive the main results of this pa- 
per: trace-norm distance upper bounds for QDD and NUDD. 
We conclude in Section [VJ and present additional technical 
details in the Appendix. 



II. MODEL 

In this section we give a formal description of the QDD and 
NUDD sequences, and the decoherence they suppress. 



A. QDD 

We describe QDD as a nesting of two UDD sequences. The 
inner UDD sequence, denoted UDD jy x , comprises N\ Z-type 
pulses, meaning instantaneous rotations by ir about the 
z-axis of the qubit Bloch sphere, i.e., Z — ie^ 1 ^/ 2 ' 1 ' 7 * . Sim- 
ilarly, the outer UDD sequence, denoted UDDjv 2 , comprises 
N2 A-type pulses, where A = ie~ 1 ^ '/ 2 ) CT * , We denote the re- 
sulting QDD sequence by QDDjVi.^- To make sure that the 
qubit state is unaltered by the sequence itself, we append an 
additional pulse at the conclusion of the sequence if N\ or N2 
is odd. Thus, if N[, denotes the number of pulses in UDD jv ! , 
where i £ {1,2}, then 



Ni = N, + Ni mod 2. 



(2) 



We define the dimensionless relative time s :— t/T, so that 
the A-type pulses are applied at times 



, • 2 3^ 

Xj = sin — — - 

3 2N ? + 2 



(3) 



where j = 1, 2, . . . , N' 2 and the Z-type pulses are applied at 
times 



A,.;, = Xj-x + {Xj - Xj-i) sin 2 



k'TT 



2N ± +2 



(4) 



where k = 1,2, . . . ,N{. 

We consider the most general form of "bounded" time- 
independent single-qubit decoherence, which is described by 
the Hamiltonian 



(5) 



where a = lg is the identity operator on the system. The 
operators B a with a 6 {0,x,y,z} are arbitrary except for 
the requirement that their sup-operator norms, i.e., the largest 
eigenvalue of (B^Ba) 1 ^ 2 , are finite 



J a ■= \\B a \\ < 00. 



(6) 



To decouple the qubit from the bath, we apply the 
QDDatj tv 2 sequence, assuming Eq. (|5}. See Fig. [T] for a 
schematic depiction of a QDD 3 3 sequence. 
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FIG. 1. (Color online) Illustration of ODD;^. Time flows from left 
to right. The switching instants are distributed according to Eqs. l|3} 
and l|4j. The pulses ir z and ir x are rotations about the z and x axes, 
respectively, by the angle ir. The total sequence time is T. 

The QDDj\Ti,jv 2 sequence is generated by the control 
Hamiltonian 

N' 2 N 2 + l K 

Hc(s) = ^(T x/ 2S(s-Xj) + ^a z ^2 X/ :s V'- W 

j=i j=x fe=i 

The corresponding control time-evolution operator is 



U c (s) = Texp C-iJ H c (s')ds'J , 



(8) 



where T is the time-ordering operator. Thus, the toggling- 
frame Hamiltonian reads 

H(s) = Ut{s)HU c {s) (9a) 
= J2 Us)o a ®B a , (9b) 

where the switching functions are 

/o(s) = l se[0,l], (10a) 

f x {s) = (-I)* 5 " 1 ae[A Jlfc _ 1 ,A Afc ) ) (10b) 

/„(«) = (-l) fc-1 (-l) J ' _1 JelAiM.Ajk), (10c) 

f z ( s ) = (-iy- 1 seix^Xj). (lod) 

For later use, we note that 

f v {s) = Ms)Us). (11) 

Next, we consider the total time evolution given by the uni- 
tary operator 

U(T)=Texp(-iT [ H(s)ds). (12) 
Jo 

Standard time dependent perturbation theory provides the fol- 
lowing Dyson series for U(T) 



00 



n=0 {a;dim(5)=n} 
-1 



(13a) 



ds n f an (s n ) / ds„_i/ a!n _ 1 (s n _i) 



rS3 r$2 

/ ds 2 fa 2 (s2) / dSif ai (si) 
JO JO 



Qs ■= (fo 



B n 



(cr ai ® B ai ) 



ra^a^a^a^ ® B a „ ■ ■ ■ B a 



(13b) 

(13c) 
(13d) 
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The notation dim(d?) = n in Eq. ( |13a| > means that the vector a 
is restricted to having n component {a n , . . . , ai}. The com- 
ponents are af. G {Q,x,y, z} for 1 < k < n. Each integer 
n a > 0, a G {0, a;, z} counts how many times the opera- 
tor <j a ® B a appears in Qg. For given dimension n of a, we 
have n — no + n x + n y + n z . In this way, the complete sum 
over all possible sequences of Bq, B x , B y , and B z are con- 
sidered. The rearrangement of the non-commuting product 
of Pauli matrices in Eq. ( |13c| > into the canonical form 
gives rise to r = (— l)* where v counts the number of times 
two different Pauli-matrices have to pass each other based on 
v a &p = ~a-pa a if a ^ (3. 

Our goal is to find an upper bound for each term fsQs 
separately. We exploit |/ Q (s)| = 1, ||er a || = 1, and 



?a|| — J a 



(14a) 
(14b) 



to obtain the upper bounding series given below for the series 
(|13a|i term by term 



oo 



\U(T)\\ < }^T n F n 

n=0 



E J * 

{a;dim(a)—n} 

Sq(T) 

exp[( J + J x + J y + J Z )T\ 



(15a) 

(15b) 
(15c) 



where 



F n ;= / ds n / ds n . 
Jo Jo 



S3 I-S2 

ds2 / dsi 
ii Jo 



1 

n! 



(16) 

The equality between ( |15a[ ) and ( |15c|> i s most easily seen by 
realizing that the right hand side of ( |15a| > is the Dyson series of 
( |15c| i as obtained by time dependent perturbation theory. This 
series will be used extensively when we compute the perfor- 



mance bounds in Section III which generalize the UDD results 
in Ref. Ifl4l. 



B. NUDD 

The NUDD sequence is a generalization of the UDD and 
QDD sequences that suppresses decoherence for a multi-qubit 
system. Assume that the system comprises m qubits. We con- 
sider the most general form of "bounded" time-independent 
?7j-qubit decoherence. Subject only to the constraint on the 
bath operators Bp that they are bounded 

Jp := \\Bgt\\< oo, (17) 

the Hamiltonian that gives rise to this general decoherence is 

(18) 



peD m 



where 



£» ro :={(0,0),(l,0),(l,l),(0,l)} m d9a) 
ap, := oi,^ ® • • • ® o m4lm . (19b) 



Here the Oj ijU . 's are Pauli matrices or the identity and we use 
j G {1, . . . , to} to index the qubits. For the jth qubit, we 
use the binary component fij 6 {(0, 0), (1, 0), (1, 1), (0, 1)} 
of the vector jl := (fii, u 2 , . . . , fi m ) to denote the Pauli matrix 
subscripts {0, x, y, z}, respectively. 

The NUDD sequence for an m-qubit system consists of 2m 
nested levels (two per qubit). Let T be the total duration of 
the NUDD sequence and N[ be the decoupling order of the 
ith-level UDD [Eq. Q]. Then the pulses at the zth level are 
applied at the instants 



\ \ • 2 hir 

- A i 2m ,b m -i,...,2 l+ Jsin 2Ar , - - 



(20) 



where 1 < i < 2m and 1 < < N[. Even values of i 
correspond to a x pulses applied to qubit number j = i/2, 
while odd values of i correspond to a z pulses applied to qubit 
number j = (i + l)/2. The control Hamiltonian is therefore 



(21) 



where the inner sums are multiple sums, one for each value 
of the index p labeling the nesting levels i,i + 1, ...,2m. The 
control time-evolution operator is 



H c {s) -- 


2m 

=- y 

2 ^ 

{i>2, even} 


K 

&i/2,x E 


S(s 




7T 

+ 2 


2m- 1 

E 

{i>l.odd} 


K 

l)/2,z E 


S(s- 





U c (s) =Texp(-i j H c (s')ds' 
Jo 



(22) 



The toggling-frame Hamiltonian is 

H(s) = U c ( S yHU c (s) (23a) 

= E ft( s ) d t® B P ( 23b > 



p£D„ 



3 = 1 



(23c) 



where the /j ijU . (s)'s are the switching functions. They can be 
obtained by the anticommutation relations of the Pauli matri- 
ces 

4(0,0) (s) = 1 (24a) 
for s e [0, 1] 

/i,(i,o)(«) = (-l)^- 1 (24b) 

forS e [k 2m ,hm-l,-,hj-l-l^hm,hm-U-,hj-l) 

/i,(0,i)(s) = (-1)^ (24c) 
for s € [M am ,h m -i,...,l2j-ii^h m ,h m -i,...,l2i) 
= /i,(i,o)(s)/j,(o,i)(s) (24d) 
fors G [0,1]. 
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In the last equation we used Eq. ( fTT| . 

Now consider the total time evolution given by the unitary 
operator 



U(T) = Texp ( - iT [ H(s)ds 
Jo 



(25) 



Standard time dependent perturbation theory provides the fol- 
lowing Dyson series for U(T) 



U(T) 



^2(-iT) n £ fp u ...,p n Qp u 

{P k ED m }% =1 



n=0 



(26a) 



fpi, 



lfp n -A s n-l) ■ ■ ■ 

(l Jo 

ds 2 fp 2 (s2) / dsif Pl (si) (26b) 
o Jo 

r J] a;*®B Pn ---B Pl , 



(26c) 
(26d) 



where r — YYjLi r j an< ^ eacn r j := (~^) Vj > wnere Vj counts 
how often two different ctj^. with the same j pass each other 
in order to obtain ( |26d| i from ( |26c| i. Each np in Eq. ( 26d| i 



indicates how many times ® Bp appears in Eq. ( |26c 
the nps satisfy J2peD m n P = n - 

Following steps analogous to the QDD case we have 



\fp\ = 1 

iqp^pj < n j ? 



(27a) 
(27b) 



peD„ 



This allows us to write 



\\u(t)\\ < j2 Tnp n E II J 7 ™ 

"=0 {p k eD m )l =1 p£D,„ 

=: Sn(T), (28b) 

where i 7 ^ was defined in Eq. ( fT6] l and the right hand side of 
d28ab is the Dyson series applied to 



(29) 



5 N (T) =exp(T J P 

P&D m 

For an alternative derivation see Section lTlIBI 



III. BOUNDS FOR GENERAL DECOHERENCE 
A. QDD 



term in the Dyson series ( |13a| i that contributes to the er a -type 
error is of order T da+1 or higher^] We will use this result in 
deriving the bounds for decoherence. 



Single-axis error a a 


Ni mod 2 


A^2 mod 2 


Decoupling order d a 


a — x 


Oor 1 


Oorl 


Ni 


a = y 




1 
1 



1 

1 


msLx{Ni,N 2 } 
max{7Vi + 1,N 2 } 
Ni 
Ni + 1 


a = z 



1 


Oor 1 
Oor 1 


iV 2 

min{iVi + 1,JV 2 } 



TABLE I. Summary of single-axis error suppression from Ref. 1181 . 

Any operator acting on the qubit subspace can be expanded 
in terms of the Pauli matrices and the identity. Hence 



U(T)= £ a a ®A a {T) 



(30) 



is another way to write Eq. ( |13a| i, which serves to define the 
bath operators A a (T). We classify the decoherence error 
based on the parities 



Pa '■= n a m °d 2, a G {x, y, z}. 



(31) 



and Using 



(T a a l3 = ie al3l <T 1 + 8 a[j t a,/3e{x,y,z} (32) 



where e Q( 3 7 is the Levi-Civita symbol, we find the results gath- 
ered in Table |lTl 



case j 


Px 


Py 


Pz 


Channel 


Decoupling order 














cro 


1 


1 








1 




d z 


2 





1 





<Jy 


dy 


3 





1 


1 


Ox 


d x 


4 


1 








Ox 


d x 


5 


1 





1 


<Jy 


dy 


6 


1 


1 





Oz 


d z 


7 


1 


1 


1 


cro 


1 



TABLE II. Classification of the error terms in the series Eq. \l3a\ 
according to the parities p a := n a mod 2. Recall that n a counts 
how many times the operator a a ® B a appears in Qs- The Channel 
column gives the corresponding single-axis error term in Eq. \2>Q) . 
The last column gives the order to which the single-axis error is sup- 
pressed. For completeness, the identity channel cro is also listed, 
though rather than being suppressed it increases in linear order. 



Similarly to the bounding series found in Ref. 11411 we can 
straightforwardly deduce bounding series for the various cases 



Very recently Ref. [ 1 8 1 found rigorous lower bounds on the 
decoupling orders of a QDDjvi,jv 2 sequence. The result is 
summarized in Table [I] We define the decoupling order of a 
single-axis error a a to be d a for a € {x,y,z}. Thus, af- 
ter applying the QDD^ jy 2 sequence, the first non- vanishing 



In Ref. 1171 the result for the decoupling order of <r z when Ni is odd was 
found numerically (using a spin-bath model) to be min{2iVi + 1, JV2}, 
which improves on the analytical bound min{A r i + 1, 7V2} from Ref. [18 ]. 
In all other cases Refs. 1 17 18] were in perfect agreement. 
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in Table |n]from Sq(T) in Eq. ( fBb) . We first write Sq(T) as 
a sum over the eight terms differing in at least one parity 



Sq(T) 



3=0 



2°p z + 2 1 p y + 2 2 Px . (33) 



The parity p a of n a is the parity of Sq(T) = exp[(J + 
Jx + Jy + Jz)T] as function of J a . Since exp(J Q T) = 
cosh( J a T) + sinh( J Q T) nicely splits even and odd contri- 
butions, we can split Sq(T) into separate bounding functions 
as listed in Table HTT1 



Px &a,x 


Py © S a , y 


Pz © $a,z 


$ ct ,p 











1 








1 








1 











1 


1 





1 











1 





1 





1 


1 








1 


1 


1 


1 



TABLE IV. The eight cases implied by Eq. ([37 



case j 


Px 


Py 


Pz 


Bounding function Sj 














e J ° T cosh(J x T) cosh(JyT) cosh( J Z T) 


1 








1 


e J ° T cosh(J x T) cosh(JyT) sinh( J Z T) 


2 





1 





e J °' r cosh(JxT) sinh( J y T) cosh( J Z T) 


3 





1 


1 


e J ° T cosh(JxT) sinh( J y T) sinh( J 2 T) 


4 


1 








e J ° T sinh( J X T) cosh( J y T) cosh( J Z T) 


5 


1 





1 


e J ° T sinh( J X T) cosh( J y T) sinh( J Z T) 


6 


1 


1 





e J ° T sinh( J X T) sinh( J y T) cosh( J Z T) 


7 


1 


1 


1 


e J °' r sinh(JxT) sinh( J y T) sinh( J Z T) 



TABLE III. Functions of the respective bounding series for the vari- 
ous contributions of parity combinations. 

We would like to use the bounding functions Sj to upper- 
bound ||j4 a ||. To this end, note first that, using Eqs. pO) and 
( |32l >, for a 6 {x, y, z} 

\ir & [<j a U(T)] = l-^2tr[(r a (Tp]Ap(T) = A a (T), (34) 



where trs denotes the partial trace over the system. On the 
other hand, using Eq. (113), 

oo 

tr s KU(T)} = J2(~iT) n Yl /s tr sKQ 5 ], (35) 

{a;dim(a)— n} 



n=0 



and 



trs[cr a Q s } = rtr[a a a^a^a^}B an ■ ■ ■ B ai . (36) 



Using the fact that only the parity pi) of the exponents mat- 
ters, we can rewrite the trace as 



*[<T a o5'Op>o2'] 

tr [ £7 .p«,ffi«o t ,«, .p»ffi«a,„ .p-ffi«a ( ,«- 



(37a) 
(37b) 



where p := (p x ,p y ,p z ) and denotes addition modulo 2. 



The values of a p are given Table IV 
Combining Eqs. (|34)-(|37), we have, 



n=0 {i3;dim(S)=n} 

where F n was defined in Eq. ([16 
In view of Table 



IV Eq. ( |38) provides the decomposition 
by parity of Eq. jl5\. Consider, e.g., the case a = x. Then 



Table IV tells us that 6 X £ 



1 only when {p x = l,p y = 

0, p z = 0} (second row) or {p x — 0,p y = l,p z — 1} (bottom 
row). In all other cases 9 x p = 0. Comparing with Table [111] 
we see that these two non-zero cases correspond to j = 4 
and j = 3, respectively. A similar argument informs us that 
9 Vt p = 1 only when {p x = 0,p y = l,p z = 0} or {p x = 

1, P y — 0,p z = 1}, which corresponds to j — 2 and j = 5 in 
and 9 z s = 1 only when {p x — 0,p y — 0,p z = 1} 
l,P y = l,p z — 0}, which corresponds to j = 1 and 



III 



Table 
or {p x 
j = 6 in Table [mj 

Now, note that the right-hand side of E q. ([38) without the 
9 a factor is just Sq(T) as defined in Eq. ( |15b| i. We can thus 
conclude from Eq. ( [38) that, due to the 9 a prefactor, its right- 
hand-side consists of only two non-vanishing terms for each 
value of 8 a , which acts like a Kronecker delta function for the 
parity triple (p x , Py ,p z ) of the bounding functions Sj, namely 
a = x =>■ j = 3, 4, a = y =>■ j = 2, 5, a = z =>■ j = 1, 6. 
Thus, 



||^(T)||<5 3 (r)+54(T) 
\\A y (T)\\<S 2 (T)+S 5 (T) 
\\A z (T)\\<S 1 (T)+S e (T). 



(39a) 
(39b) 
(39c) 



Next, to account for the suppression of decoherence by 
QDD we define 



Pk 



fe! dT k 3 



T=0 



so that 



Si 



(40) 



(41) 



k=0 



We consider the partial Taylor series A?' of the analytic 
bounding functions Sj which leave out the contributions up 
to and including T d 



A 



OO 

E pf Tk 

k=d+l 



(42) 



According to Tables [II] and III the contributions of case j 
are bounded by the corresponding partial Taylor series 
where d is the decoupling order of the channel. This is a man- 
ifestation of the fact that the QDD sequence causes the first d a 
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powers in T of A a to vanish, i.e., A a (T) = <D{T d « +1 ) lfl8l . 
Using Eq. ( [39) 1 we deduce that 

L X :=A^+A^>\\A X (T)\\ (43a) 
L y :=A^+A^>\\A y (T)\\ (43b) 
:= A^+A^ > P 2 (T)||. (43c) 

This is our first key result for QDD. 

Due to the analyticity in the variable T of Si for each j, we 
know that the residual term vanishes for d — > oo, that is 

lim A { [> = 0. (44) 

d— ¥oc> 

We define the dimensionless parameters (h is set to unity) 

e := J T Va := J a /J (45) 

where a € {x,y,z}. Thus, we can write the bound- 
ing functions Sj as Sj(e,ff), where ff = (r] x ,T] y ,r) z ). By 
Taylor-expanding Sj(e, r) X) r) y ,r) z ) we can express the bound- 
ing functions in terms of e and functions {g\ \ (ff)}n=i given 
in Appendix [A] 

oo 

Ai £= E 9n\n)e n (46a) 

n=d a + l 

= g d 3) a+ M^ da+1 +0(e d - +2 ). (46b) 

The behaviors of the bounding functions L x , L y and L z are 
shown in Fig. [2] for different ff. Here we assume that the de- 
coherence is isotropic, i.e., 77 := rj x = r) y = r\ z . Given this 
assumption, it makes sense to choose N\ = N 2 so that the 
decoupling orders for different types of errors are close. To 
simplify the computation, we choose both Ni and N2 to be 
even and thus, according to Table [I] d x = d y = d z . Note that 
in this situation L x , L y and L z are identical. 
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FIG. 2. (Color online) The upper bound L a for ||^4:c|| = \\A y \\ = 
\\A Z \\ according to Eqs. (|43t in the isotropic case, i.e., r\ x = r\ y = r\ z , 
is given as a function of e for various number of pulses Ni = N% = 
N G {2, 6, 16, 34}, at fixed values of r/ a . In each panel the curves 
become steeper as iV increases. 



We also investigated the case where rj x = rj y 7^ r\ z and how 
the parities of N\ and N2 affect the bounds. In Fig. [3] and 
Fig. [4] we plotted two cases. The first is where both Nx and 
N2 are even, which are represented by the thick lines. The 
second is where Ni and A^2 are either even or odd, which are 
represented by the thin lines. We picked r\ z to be 10 -2 for all 
the plots and rj x = t] y G {10~ 4 , 10~ 2 , 1, 10 2 }. Since r\ z is 
fixed, we also fixed the values of N2 for the two cases. We 
picked N2 to be 10 in the first case and 9 in the second case to 
explore the effect of parity. We varied the value of Ni based 
on the values r/ x and r\ y . We can see from the figures that in 
the majority of cases, the parities of N\ and N2 only change 
the decoupling order by 1 or do not change it at all, so the 
bounds do not vary much. However, when N\ is small com- 
pared to N 2 , the parities can change the bounds nontrivially 
as seen in the top left panels of Fig. [3] and Fig. [4] We can 
also see this from Table[I] When A^i is small compared to N2, 
the decoupling orders of a y and a z can decrease considerably 
when Ni switches between even and odd values. 
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FIG. 3. (Color online) The upper bounds for ^A x ^ (dotted lines), 
\\A y \\ (dashed lines) and \\A Z \\ (filled lines) according to Eq. l |43| l 
are depicted. Here ^2 is 10. The values of iVi in the four panels 
going from left to right and top to bottom are 2, 10, 18 and 34. The 
value of r\ z = 10 -2 . 

Equation ( |43] ) allows us to establish tight bounds on the rep- 
resentation of the unitary time evolution as given in Eq. (30) . 
In subsection llV Al we use this fact to establish a bound on the 
trace-norm distance between the unperturbed qubit state and 
the protected qubit state. 

B. NUDD 

If we were to use the same method for NUDD for m qubits 
as we did in the QDD case for one qubit, we would need 
to consider 2^ 3 ) cases based on the parities of np. Re- 
call that np indicates how many times ap Bp appears in 
ftiln • • • (ap 1 ^iBpA; Eq. ( |26d) . To simplify the anal- 

ysis, unlike the QDD case, we do not address each error term 
separately. 

Any operator acting on the space of m qubits can be 
expanded in the {<rp} basis. Generalizing the QDD case 
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4. (Color online) The upper bounds for || Aa; || (dotted lines), 
alj (dashed lines) and (filled lines) according to Eqs. {43} 



are depicted. We make A?2 to be 9 in this case. The values of N\ in 
the four panels going from left to right and top to bottom are 3, 10, 19 
and 34. The value of r\ z — 10~ 2 as in the previous case. 



[Eq. ( [30] )] we can thus write 

u(t)= y, $p®Mn 



(47) 



p£D„ 



where ap and D m were defined in Eqs. (fT 9l. Eq . ( |47] i provides 
another way to organize the terms in Eq. ( |26a| i. We introduce 
the upper bounding series Sn for each Ap such that 

S N (T) = Y S p( T ) 
peD n , 



(48a) 
(48b) 



holds. Below we construct this series explicitly. Let 

0:={0,...,0} (49a) 

K := D m \0. (49b) 

Since is the all-identity term, it is the term that causes no 
errors; its corresponding bounding series is Sg(T). On the 
other hand, every ap with jl € K, is an error term with corre- 
sponding bounding series Sp(T). For later use, we note that 
S g (0) = 1 and S p (0) = for all fx G K. 

In order to derive bounds we first recall that U (T) from 
Eq. ( p5| ) is the solution of the Schrodinger equation 

d T U{T) = -iH{T)U{T) (50a) 



]T fp{t/T)dpBp 

p£D m 



Its integration from Sn(0) = 1 recovers the previous result 
( |29"1 ) precisely. 

In analogy, the differential equation for each ^(T) reads 



d T Ap(T) 



fp®v(t/T)Rp(^p-pBp^c; 



u£D„ 



Av{T), 



(52) 

where, as before, © stands for sums modulo 2. Note 
that addition and subtraction are equivalent modulo 2: fl © 
V = jl Q v. It is essential that the sums modulo 2 in 
{(0, 0), (1, 0), (1, 1), (0, 1)} faithfully reproduce the spin al- 
gebra of the identity and the Pauli matrices except for factors 
of ±i, which are collected in Rp^.^. Replacing operators 
by their norms and complex numbers by their moduli on the 
right hand side yields the generating differential equation for 
the bounding series Sp satisfying Eq. d48ab, 



d T Sp{T) 



PeD n 



Sa(T). 



(53) 



Unfortunately, the set of differential equations ( |53) is still 
too difficult to be solved generally. Hence we aim at a looser 
bound by defining 



Jl 



:= max Jp 
p&ic 

= J l 



(54a) 

(54b) 
(54c) 



If we substitute J\ for each Jp with jl € JC we obtain one 
bounding series S X {T) for all Sp(T) = S X {T). To see this 



U(T) (50b) 



with U(0) = 1. Eq. ( |50b| i generates exactly all terms of the 
series of U (T). Thus, in order to obtain the series for S^(T) 

all we need to do is to replace — i ->• | — i\ = 1, fp -> |/^| = coiTesponding matrix are A + 
1, dp ->• \\a~p\\ = 1, and Bp ->• \\Bp\\ = Jp on the right 
hand side. This provides us with the generating differential 
equation for the bounding series 



we insert Si(T) = Sp(T) on the right hand side of Eq. ( |53| ) 
yielding 

d T S 5 (T) = J g S 6 (T) + yJiSi(T) (55a) 
d T Sp(T) = J ff Si(T) + JiS s (T) + ( 7 - 1) JxSxiT) (55b) 

where jl € JC and 7 := 4 m — 1 = \JC\ is the number of 
non-identity terms. The first term on the right hand side of 
( |55a| i results from = 000, the second from v = © V for 
v £ JC. The first term on the right hand side of ( |55b| i results 
from = jl © jl, the second from jl = ft © 0, and the third 
from jl® V where v £ JC\fl. There are 7 — 1 terms of the 
latter kind, independent of jl. This independence implies the 
equality of all Sp(T) for jl £ JC because they all start at the 
same value Sp(0) — 0. Thus, replacing Eq. ( |55b| i, we have 

drSi(T) = J S Si(T) + Ji%(T) + (7 - l)JiSx(T) (56) 



p£D m 



which, together with Eq. ( |55a[ ), constitutes a two-dimensional 
set of linear differential equations. The two eigenvalues of the 

r o + "fJi and A_ = Jq - Ji 
so that the general solution reads S\{T) := Aexp(X + T) + 
B exp(A_T) with coefficients A and B. The initial conditions 
Sg(0) = 1 and 5 X (0) = imply 9 T 5i(0) = J x from which 
^4 and B are determined to yield 



5 N (T). 



(51) 



5i(T) 



exp(J g r) 
7 + 1 



exp( 7 J 1 T) - exp(- JiT)). (57) 
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The bounding series of the sum of all error terms is given by 
Sk;(T) := jSi(T) which reads 



Sk(T) = 7eXP [^ T) (exp( 7 J 1 T) - exp(-J 1 T)). (58) 
7 + 1 



This result is the bounding series for all error terms in NUDD. 
We define 



Pk 



1 d k 

OO 

E 

k=N+l 



Sk 



PkT 



T=0 
k 



(59a) 
(59b) 
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FIG. 5. (Color online) The upper bound Ad mln for X^gjc ll^ll 
is given as a function of e for various pulse numbers d m i n G 
{5, 10, 20, 40}, at fixed values of rj, and for m — 10. In each panel 
the curves become steeper as d m i n increases. 



d min := min {N' t } (60) 

l<K2m 

be the decoupling order of the NUDD sequence |[T9l . Since 

H^rll < %(T),wehave 

E H^ll < E %( T ) ^ S ^ T )> (61) 

and hence 

]T H^ll <A dmin . (62) 

By Taylor-expanding Eq. ( |58| we can write Sjc (T) and thus 
Ad as an infinite series in terms of e and r\ of Eqs. ( |54c| i and 
(54b]). We do not write the result explicitly, but the formula 
(58|l for Skz(T) is simple enough to allow the bounds to be 
computed easily with any computer algebra program. This is 
our first key result for the NUDD sequence. 

To leading order in e, we find 

A d = g d+1 (r 1 ,m)e d+1 +(D(e d+2 ), (63) 

where 

,h, m ):=(l- W 1 - + «-;>'-' 1 -< (64) 

In Fig.[5]we include plots of Ad mto as a function of £ for var- 
ious values of <i min and rj. We can see that, as expected, the 
behavior of Ad min is similar to that of L a in the QDD case. 



IV. DISTANCE BOUND 
A. QDD 

Next, we shall use the bounds on the bath operators A x , A y 
and A z to derive a bound on the trace-norm distance 

D[p s (T),p° s (T)] = l\\p s (T)- pKT)^ (65) 

between the actual qubit state 

ps(T) :=tr B [p SB (T)], (66) 
and the "error-free" qubit state 

^(T):=tr B [p° B (T)] (67) 

where p^ B (T) is the time-evolved state of the whole qubit- 
bath system without coupling between the qubit and the bath, 
and trg denotes the partial trace over the bath degrees of free- 
dom. The trace norm \\A\\i is the trace of (A^A) 1 / 2 . The 
trace-norm distance D is the standard distance measure be- 
tween density matrices (TJ. The method we employ here is 
similar to the one in Ref. [ 14 1. We consider an initial state 

pis(0) = hWI®PB (68) 

in which the qubit is in a pure state and the bath is in an 
arbitrary state p B (e.g., a mixed thermal equilibrium state). 
The initial state evolves to psb(T) = C/(T)pg B (0)C/ t (T) 
when the qubit and the bath are coupled, or to p$ B (T) — 

[l s ® U B (T)]p% B (0)[t s ® Ul(T)] when the qubit is isolated 
from its environment. The unitary bath time-evolution opera- 
tor without coupling reads 

U B (T) := exp(-iTB ) (69) 

where Bq is the bath term in Eq. Q. 
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Let us first define the bath correlation functions 

b a0 (T) :=tr[A a (T)p B AUT)} 



(70) 



where a, ft € {0, x,y, z}. To simplify the notation, we will 
omit the time dependence T in b a p (T) and A a (T) from now 
on. Computation yields (see Appendix |B \\ 



+\ E (IM-M). (7i) 

We know from the unitarity of U and Eq. d30l) that 



1 = U^{T)U{T) 
= l s <gi E A^ a A a 

a£{Q,x,y,z} 

+ E ° a ®{AiA + AlA a + i[Al,A 7 }), 

a£{x,y,z} 

(72) 

where in the last sum 7} are adjusted to a so that {a, f3, 7} 
is a cyclic permutation of {x, y, z}. Therefore we have 



1 


= AlA - 


h 4^ - 


f Ay Ay - 


f- a\a z , 


(73a) 





= 4a>- 


h 4^ - 


*- *4 A * 


- iA z A y , 


(73b) 





= A\A Q - 


1- A\Ay - 


h i4^ 


— iA x A z , 


(73c) 





= A\A - 


h a\a z - 


- iA x A y 


— iAyA x . 


(73d) 



D[ Ps (T),p° s (T)} 
< \\AJ\ + \\A. 



+ HAr|||K|| 



A 



\A Z 
I A, 



(2) 



\\A Z 



hi 



\Ayf 

-\\\\A*\\ 

L X Ly 



\A Z \\ 2 
(77a) 



LyL z 



A 



(5) 



(2h2 



,(6) 



(77b) 
(77c) 

,(!h2 



A^Aj) 



AipAi 5) + A?>A 



(4) A (2) _ 



A^Ag) 
A^Ag) 



a£>a£> 
a^a« 



H2A«A^ 
A^Aj 

A^Ag. 



This is the rigorous bound on the trace-norm distance and 
hence our second key result for QDD. 

Using Eq. ( |46b| i, one can rewrite this bound in terms of e 
and fj, which we do not write out here for brevity. 

In explicit calculations the leading order in e will dominate 
for e — > 0. Then only the first line in Eq. ( |77c| > plays a role, 
since all other terms are of higher order. Hence we have 

D[p s (T), P ° s (T)} < [gl } + M+9^ +i m^ +1 



Equation ( |78] l is our third key QDD result. It shows that the 
QDD bound is dominated by the channel with the smallest 
decoupling order d a , as was expected intuitively. 



It follows from Eq. ^ that (i\ £ q6{0 i ^ z) A^A a \i) = 1 
for all normalized states \i), and thus 

(t\AiA \i) = \\A \ l )\\ 2 <l (74) 

because (i\ T, ae {x,y,z} A l A a \i) = Eae{*,w,*} \\ A a\i)f is 
nonnegative. In particular, we know that 

Poll < I- (75) 

To obtain a bound on the functions 6 Q( g we use the following 
general correlation functions inequality (see Ref. lfl4l for a 
proof) 

HQp B Q'}\ < WQWWQ'l (76) 

which holds for arbitrary bounded bath operators Q,Q'. Us- 
ing Eq. (|75j , Eq. f75) , and the bounds (|45]l in Eq. {71} yields 



B. NUDD 

We consider an initial state 

p° SB (0) = \i>)(ij\®PB (79) 

where \ip) is the state of the m-qubit system. The ini- 
tial state evolves to p SB {T) = U (T) p" SB (0)W (T) when 
the qubit and the bath are coupled, or to p^ B (T) — 
[Is ® Ub(T)]p% b (0)[± s ® Ul(T)] when the qubit is isolated 
from its environment. The unitary bath time-evolution opera- 
tor without coupling reads 

17 B (T) := exp(-iTB s ) (80) 

where is the bath operator that corresponds to the identity 
system operator. We know from the unitarity of U (T) and 
Eq. ((47]) that 

E 44 = 1- (81) 
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Following the same steps as in the QDD case yields 



\A 6 \\ < 1. 



(82) 



Explicit computation (see Appendix B 2 1 as in the QDD case 
shows that 

D[p s (T),p° s (T)} < ||^||) 2 + Poll £ WM • < 83 > 

We use the bounds on the operators Ap to derive a bound on 
the trace norm distance between psb{T) and p$ B (T). Substi- 
tuting Eq. (|62|) and Eq. ([82} into Eq. (|83|l, we conclude that 



(84) 



D[p s (T),p°(T)]<A2 +A dmin , 



where Ad min is the bound presented in Eq. ( |59b| i. The bound 
is dominated by A<j min which is the leading order term for 
e 0. By expanding Ad min we find for the leading order 
term 

D[p s (T),p° s (T)} < g dma+1 ( V ,m)e d ™ +1 +0(e d ^ +2 ). 

(85) 

Equation (f85]> is our second key result for NUDD. 



V. CONCLUSIONS 

We have derived and presented rigorous performance 
bounds for the QDD and NUDD sequences, respectively. 
These sequences, which build on the UDD sequence, protect 
a qubit or system of qubits from general decoherence, under 
the assumptions that the pulses are instantaneous and the bath 
operators are bounded in operator norm. Our key bounds are 
given in Eq. {77c) for QDD and in Eq. <[84j» for NUDD. The 
leading order terms are identified in Eq. (|78[) (for QDD) and 
Eq. ((83]) (for NUDD). These results show that if the total se- 
quence time is fixed, we can make the error D[ps(T), p^(T)] 
arbitrarily small by increasing the number of pulses at each 
UDD level comprising the QDD or NUDD sequences. 

When instead the minimum pulse interval is fixed, we ex- 
pect that, just as in the case of UDD lfT4l . there will be an 
optimal sequence order, beyond which performance starts to 
decrease. A rigorous study of this aspect of QDD and NUDD 
is an interesting topic for a future investigation. We hope that 
the results presented here will inspire experimental tests of 
the QDD and NUDD sequences in physical systems with bath 
spectral densities exhibiting relatively hard high frequency 
cutoffs, a condition which corresponds to our key assumptions 
[Eqs. (j6j,([T7j] of bounded bath operators. 
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Appendix A: Bounding polynomials 

The polynomials that appear in the Taylor expansions of the 
bounding functions Sj(e,ff) are: 



g^(ff) = 




F Vx ~ 




+ v z ) 


+ (i- 


- Vx - 


VVy- 






+ (H 


- Vx - 


Vy - 


\ 1 


+ (i- 


Vx - 


Vy + 


vzY 




-(14 


- Vx + 


Vy - 


-vz) 1 


-(i- 


Vx + 


Vy- 


vzY 




-(14 


- Vx - 


Vy ~ 


-vz) 1 


-(i- 


Vx - 


Vy - 


vzY] 

(Al) 


g?\v) = 


srK 1 


F Vx - 


IVy 


+ %) 


+ (i- 


- Vx - 


VVy - 


^vzY 




-(14 


' Vx - 


Vy " 




-(i- 


Vx - 


Vy + 


VzY 




+ (M 


- Vx + 


Vy - 




+ (i- 


Vx + 


Vy - 


vzY 




-(M 


- Vx - 


Vy - 


-vz) 1 


-(i- 


Vx - 


Vy - 


vzY] 

(A2) 


gl 3) (v) = 


1 


F Vx - 


VVy 




+ (i- 


- Vx - 


VVy - 


^vzY 




-(14 


' Vx - 


Vy~ 


^vJ 


-(i- 


Vx - 


Vy + 


vzY 




-(14 


- Vx + 


Vy - 


-vz) 1 


-(i- 


Vx + 


Vy - 


vzY 




4-(l4 


- Vx - 


Vy " 


-v z y 


+ (i- 


Vx - 


Vy ~ 


vzY] 

(A3) 


rf 4) 0) = 


^ 


F Vx - 


VVy 


+ %) 


-(i- 


- Vx - 


\-r)y - 


^vzY 




+ (14 


- Vx ~ 


Vy - 


^v*Y 


-(i- 


Vx - 


Vy + 


vzY 




+ (14 


- Vx + 


Vy - 


~vzY 


-(i- 


Vx + 


Vy - 


vzY 




+ (14 


- Vx - 


Vy - 


-vzY 


-(i- 


Vx - 


Vy ~ 


vzY] 

(A4) 


gl 5) (v) = 


8^ 


F Vx - 


VVy+Vz) 


-(i- 


- Vx - 


VVy - 


^vzY 




+ (14 


- Vx - 


Vy ~ 


*-vzY 


-(i- 


Vx - 


Vy + 


vzY 




-(14 


' Vx + 


Vy - 


-vzY 


+ (i- 


Vx F 


Vy ~ 


vzY 




-(14 


- Vx - 


Vy " 


-vzY 


+ (i- 


Vx - 


Vy ~ 


vzY) 



(A5) 



11 



9l e Hff) = -[(l + Vx - 


'rVy + Vz) 


-(1- 


" Vx - 


\-Vy 


1- Vz 


-{l + Vx- 


Vy + VzY 


+ (1- 


Vx - 


Vy + 


Vz) 


+ (l+Vx + 


Vy - Vzf 


-(1- 


Vx + 


Vy ~ 


Vz) 


- (1 + Vx - 


Vy - VzY 


+ (1- 


Vx - 


Vy ~ 


Vz) 



I 
I 

il 

(A6) 

Appendix B: Distance Bound Calculation 
1. QDD 

In this appendix we prove the bound on the trace-norm dis- 
tance in Eq. (jTTJ. To simplify the notation, we let 



Pn 



(Bl) 



Then 



2D[p s (T),p° s (T)} 
= \\tr B [p SB (T)-p° SB (T)] || 1 

= |itb [U(T)p° SB (0)tf(T)] -tr B [p° SB ( 



tr B 



E ct ° 

aG{0,a;,y,2} 



(Po®Pb) 



a£{0,a;,?/,2:} 

|| (&oo - l)-Po + b 0x P a x + b 0y P a y 
+ bo z Pocr z + b X Q<j x Po + b xx a x Poa x 
+ b xy a x Po<jy + b xz a x Po<j z + b y o<j y Po 
+ b yx a y Po<r x + b yy a y Po<r y + b yz a y Poa z 
+ b z o<J z Po + b zx a z P a x + b zy a z P a y 
+ b zz a z P a z \\ 1 

-\b Qy \ + \b Qz \ 
\b xy \ + \b xz \ 
\b yy \ + \b yz \ 
\b zy \ + \b zz \. 



< l&DO — 11 + M 
+ |&sco| + \bxx\ ■ 

+ IM + \by x \ - 

+ \b z0 \ + \b zx \ - 



(B2a) 



(B2b) 



In going from Eq. ( |B2a| > to Inq. ( |B2b| i, we used the triangle 
inequality, the unitary invariance of the trace norm, and the 
normalization of This proves Eq. ( f7l) , 
We require one last result: 



|6oo-l| = |tr[^ p B ^o] 



1| 



HaIa qPb ] - 1| 

|tr{[l - A\A X - A y A y - A\A z ]p B } - 1| (B3) 
|tr[pe] - tr[A xPB Ai] - tr[A y p B Al] - tr[A z p B A\}\ 
\b xx ~t~ b yy -\- b zz \, 



where we used cyclic invariance of the trace in b a p together 
with Eq. ( |73a| i and the normalization tr[p B ] = 1. Eq. ( fTTj i now 
follows immediately from the triangle inequality. 

2. NUDD 



Here we prove Eq. < |83p following the same steps as in the 
QDD case. 



D[p s (T),p° s (T)} 
= \\Mpsb(T)-p° sb (T)}\\ 1 

= \\MU(T)p° SB (0)u^T)} - HbMo)]^ 



tr B 



E a P A P 



(Po®Pb) 



E «?4 

P&{0,1} 2 ™ 



Po 



E KAlpsA^aflPoav - P 

P&{0,l} 2m £?6{0,l} 2m 

< |tr(^p B ^) - 1| + J2 £ M A ]tPB A *)\ 
+ £ MAlpBAp)\ + J2 H4pbA 6 )\ 

< £ itr(4p B ^)i + E E m a Ipm\ 

+ J2 HAIpbA p )\ + J2 |tr(4p B ^)| 

< ( E ii^ii) 2 + E h^ii 2 + 2iiA 5 n \\M ( B4a > 



p£lC 



peic 



pete 



< 



2 (EH^Il) 2 +2||A H|^||^| 



(B4b) 



p£K 



peic 



In the derivation we used the unitary invariance of the trace 
norm, triangle inequality, the normalization of Eq. ( |76*| ) 
and Eq. ( |8Tj l. The step from inequality ( |B4a| > to ( |B4b) > is not 
necessary, but loosens the bound for the sake of simplicity of 
the result. 
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